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Abstract 

We present the relation between the atomic and spin-electronic structures of infinite single-row 
atomic wires made of Al atoms during their elongation using first-principles molecular-dynamics 
simulations. Our study reveals that the Peierls transition indeed occurs in the wire with magnetic 
ordering: it ruptures to form a trimerized structure with antiferromagnetic ordering and changes 
from a conductor to an insulator just before forming a linear wire of equally-spaced atoms. The 
formation of the trimerized wire is discussed in terms of the behavior of the u-symmetry bands of 
the Al wire. 

PACS numbers: 73.21.Hb, 75.75.+a, 68.65.La, 73.63.Nm 
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In the past decade, there has been a great deal of interest in the preparation of atomic 
wires and the investigation of their mechanical and electrical behaviors |IJ. The properties 
of such minute systems with low dimensionality can be different from those of the bulk. 
In particular, an infinite single-row atomic wire (ISAW) has attracted great attention from 
a fundamental point of view, because phenomena peculiar to the reduced dimensionality 
are expected to emerge best for the ISAW having an ultimately simplified one-dimensional 
structure |2], ||, £|, ||, ||. Half of a century ago, Frohlich and Peierls |§ discovered instability 
in one-dimensional electron systems that has come to be called the Peierls instability: the 
regular chain structure in a one-dimensional wire with a partly filled band will never be 
stable, since one can always find a more suitable structure, such as a dimerized structure, 
for which a break will occur at or near the edge of the Fermi distribution. Thus the band 
will be split into a number of smaller zones, some of which are filled and the rest empty, and 
as a consequence, the wire changes from a conductor to an insulator, which is the so-called 
Peierls transition. So far, several quasi-one- dimensional metals such as organic ones were 
observed to exhibit the Peierls transition ||. 

First-principles molecular-dynamics (FPMD) simulations have been employed to eluci- 
date the atomic structures and spin-electronic properties of atomic wires || |], ^ The 
studies for an ISAW made of gold atoms ||, |G| demonstrated that the Peierls transition 
certainly occurs: the gold wire breaks to form a dimerized structure accompanied with a 
metal-insulator transition when the average interatomic distance increases beyond ~ 5.7 
a.u., where a zigzag geometry transforms into a linear one in the process of its elongation. 
What has to be noted is that a gold atom is considered simple due to its s-type outermost- 
valence orbital and its nonmagnetic property. The situation may be more complex in the 
case of monatomic wires made of elements having valence orbitals other than the s-type one 
and/or made of magnetic elements. Recently, for an ISAW consisting of aluminum atoms 
(Al-ISAW) with 3s- and 3p-type valence orbitals, Ayuela et al. || reported that an ISAW 
with a regular chain structure of equally-spaced atoms on a straight line is in a metallic 
state with certain magnetic ordering and it is energetically more favorable than a dimerized 
wire at an interatomic distance of 5.3 a.u., which is close to the nearest-neighbor atomic 
distance of 5.4 a.u. in the aluminum bulk. Thus, the existence of a metallic linear ISAW of 
equally-spaced atoms gives rise to a puzzle in Peierls's view that such an atomic wire will 
never be stable, although the result is appealing in that a stable ISAW is in a magnetic state. 
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Another theoretical study by Sen et al. Q suggested the possibility of the Al-ISAW to rup- 
ture into a tetramized wire under the assumption that the effect of the distortion on cohesive 
energy is negligibly small, but they overlooked the occurrence of magnetic ordering in the 
wire. Therefore it is imperative to clarify the properties of the Al-ISAW for understanding 
of the Peierls instability in a one- dimensional atomic wire during its elongation. 

In this Letter, we implement detailed FPMD simulations to elucidate the close relation 
between the atomic and spin-electronic structures of the Al-ISAW. Our finding is that the 
Peierls transition indeed occurs in the Al-ISAW with magnetic ordering: the ISAW during its 
elongation breaks just before it transforms from a zigzag structure to a linear one. When the 
wire fractures, the spin-density wave (SDW) arising from the Coulomb repulsion interaction 
leads the wire to spontaneously stabilize, and opens up an energy band gap at the Fermi level 
(Ep). Consequently, a trimerized wire with antiferromagnetic (AFM) ordering emerges, i.e., 
the wire transforms from a conductor into an insulator at the formation of the trimerized 
wire with AFM ordering. Moreover, we confirm that this distorted structure is dominated 
by the behavior of the cr-symmetry bands crossing the Ep. 

Our FPMD simulation is based on the real-space finite-difference method [H|, which en- 
ables us to determine the self-consistent electronic ground state and the optimized atomic 
geometry with a high degree of accuracy by making use of the timesaving double-grid tech- 
nique ]IT| and the direct minimization of the energy functional JP2 |. The norm- conserving 
pseudopotentials |T3| of Troullier and Martins [[H are employed to describe the electron-ion 
interaction. Exchange-correlation effects are treated with the generalized-gradient approx- 



imation |15j in the density-functional theory. The nine-point finite-difference formula (i.e., 
the formula with iV=4) is adopted for the derivative arising from the kinetic-energy 
operator. We take a cutoff energy of 27 Ry, which corresponds to a grid spacing of 0.6 a.u., 
and a higher cutoff energy of 247 Ry in the vicinity of the nuclei with the augmentation of 



double-grid points |TT] . We employ a technique that involves the use of a periodic supercell 
whose size is chosen as L x =L y =20 a.u. and L z =N atom x d av . Here, the z axis is taken along 
the wire axis, the x,y axes are perpendicular to the wire, L x , L y and L z are the lengths 
of the supercell in x, y and z directions, respectively, N atom is the number of aluminum 
atoms within the supercell, and d av is the average value for projections of the interatomic 
distances between adjacent atoms onto the z component. In order to eliminate completely 
unfavorable effects of atoms in neighbor cells which are artificially repeated in the case of 
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FIG. 1: (A) Bond angle 9 (closed circles) and interatomic distances d\ (crosses), (B) restoring 
force (closed circles) and total energy per atom (crosses), (C) magnetic moment for the model of 
the Al-ISAW as shown in the inset of (A). The supercell of N a t om =2 is employed. The zero of 
energy is chosen to be the total energy of an isolated aluminum atom. In the inset of (A), the 
closed circles denote the atoms in a supercell and the open circle shows the replicated atom in the 
adjacent supercell. 

periodic boundary condition, we impose the nonperiodic boundary condition of vanishing 
wave function out of the supercell in the x and y directions and adopt the periodic boundary 
condition in the z direction. 

First we examine the stability of the Al-ISAW in the case of N atom =2. The calculation 
model, the parameters of the fully optimized geometry of the ISAW, the restoring force, the 
total energy per atom, and the magnetic moments are depicted as a function of d av in Fig. [I], 
and the energy band structures of the ISAW are plotted in Fig. Q In these calculations, 24 
k z points are employed for the integration over Brillouin zone. As an initial configuration 
for FPMD calculation, the parameters of the atomic geometry d\ and 9 at d av =4A5 a.u. are 
randomly set in the ranges of 4.45±0.05 a.u. and 140±5 degrees, respectively E|. Then, 



4 



FIG. 2: Energy band structures of the Al-ISAW at (A) d av =5.2 a.u. and (B) d av =5.9 a.u. plotted 
along the direction k z in reciprocal space. The supercell of N atom =2 is employed. The zero of 
energy is chosen to the Fermi level Ep- Solid and dashed curves represent up-spin and down-spin 
electron bands, respectively. The 7r-symmetry bands are doubly degenerate. 

we increase d av and relieve the force on atoms repetitiously. In each case, the structural 
optimization is performed until the remaining forces acting on the atom are smaller than 
82.4 pN. From Figs. |I| (A) and ^|, one sees that a zigzag wire changes to a linear one upon 
elongating d av more than ~ 5.2 a.u., and even when further stretching up to d av =h.Q a.u., 
the wire does not form a dimer structure, but keeps a regular chain structure of equally- 



spaced atoms on a straight line [17]], and it remains a conductor. These results throw us into 
a puzzling in the view of the Peierls instability. The wire in ferromagnetic (FM) state has 
been also observed to be energetically more stable than those in other magnetic states, i.e., 
AFM and paramagnetic (PM) orderings, over the range of d av calculated here in the case of 

N atom = 2 • 

To gain an insight into which structure is the most stable after the Al-ISAW breaks, we 
evaluate the total energies of the wires with several atomic and spin-electronic structures. 
Whereas the atoms in the wire were alowed only to dimerize in the above study, additional 
geometrical degrees of freedom, such as trimerization and tetramerization, are now taken into 
account for definite clarification of the atomic and spin-electronic structures of the Al-ISAW 
during its elongation. The calculation models and results are depicted in Fig. ||. We take 
the average interatomic distance d av of 5.6 a.u., which is longer than the nearest-neighbor 
atomic distance of 5.4 a.u. in the fee aluminum crystal and than the critical distance of 5.2 
a.u. where an ISAW transforms from a zigzag structure into a linear one. The Brillouin-zone 
integrations are performed using 12, 8, and 6 k z points for the models (A), (B), and (C) with 
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FIG. 3: Models for (A) dimerized wire of N a t om = 4, (B) trimerized wire of N a t om = 6, and (C) 
tetramized wire of N a t om = 8. (D) Total energy per atom as a function of the bond length of the 
dimer in the wire with PM, FM, and AFM orderings. The zero of energy is taken to be the total 
energy of the wire with PM ordering at Ad=0 a.u. The average interatomic distance d av is set to 
be 5.6 a.u. In (A), (B), and (C), the closed circles denote the atoms in a supercell and the open 
circle shows the replicated atom in the adjacent supercell. 



supercells of N atom =4:, 6, and 8, respectively. In Fig. | (B), the middle atoms within the 
trimers are located at the centers of the respective trimers in an intertrimer distance of 3d av . 
In the case of the tetramized ISAW in Fig. ^ (C), the two middle atoms inside the tetramers 
are relaxed through the FPMD simulation, while both the edge atoms of the tetramers are 
frozen so as to maintain a tetramer length of 3(d av — Ad) and an intertetramer distance of 
4d av . The electronic ground states with PM, FM and AFM orderings are observed, but that 
with the ferrimagnetic ordering is not found. In contrast to the conclusion drawn by Sen et 
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FIG. 4: Energy band structures of the Al-ISAW: (A) dimerized wire at Ad=0.0 a.u. with AFM 
ordering, (B) trimerized wire at Ad=0.5 a.u. with AFM ordering, (C) tetramized wire at A<i=0.6 
a.u. with AFM ordering, (D) dimerized wire at Ad=0.0 a.u. with FM ordering, (E) trimerized 
wire at Ad=0.5 a.u. with FM ordering, and (F) tetramized wire at Ad=0.4 a.u. with FM ordering 
plotted along the direction k z in reciprocal space. Same as Fig. |3]for the definition of Ad. In (D), 
(E), and (F), solid and dashed curves represent up-spin and down-spin electron bands, respectively. 
The zero of energy is taken to the Fermi level Ep. The average interatomic distance d av is set to 
be 5.6 a.u. 

al. we can find a sensible difference in energy between regular and distorted wires. The 
total energy per atom of the trimerized wire with AFM ordering is the lowest, as shown in 
Fig. |3| (D). Therefore we conclude that the trimerized Al-ISAW in AFM state is the most 
likely structure for the distorted Al-ISAW. 

We next consider why the Al-ISAW prefers a trimerized structure with AFM ordering 
after it breaks. This reason is understood by noting the behavior of the cr-symmetry band. 
In Fig. 0, one see that the stretched wire has doubly-degenerate 7r-symmetry bands and non- 
degenerate cr-symmetry ones near the Fermi level Ep, and the 7r-symmetry and cr-symmetry 
bands cross the E F near k z w ^— and k z w ^— , respectively. When the IS AW ruptures, 
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TABLE I: Total energy per atom (eV/atom) of (A) the trimerized Al-ISAW with AFM ordering 
and (B) the Al-ISAW calculated using the model in Fig. |l| (A). The zero of energy is chosen to be 
the total energy of an isolated aluminum atom. 
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the SDW which has Qd av periodicity appears so as to open up an energy gap of cr-symmetry 
bands at the Ep and to lower the total energy. Figure |] (B) shows the energy band struc- 
tures of the trimerized ISAW with AFM ordering at d av =5.6 a.u., for which one can find 
the energy band gap at the Ep. On the other hand, we can confirm that the SDWs which 
have the other periodicity, such as 2d av and 3d av , lead the cr-symmetry bands to cross the 
Ep and increase the total energy of the ISAW relative to the trimerized ISAW with AFM 
ordering having SDW of 6d av periodicity (see Fig. [|). 

Finally, we explore the critical interatomic distance where the Al-ISAW transforms from 
the zigzag structure to the trimerized one during its elongation. Table | shows the total en- 
ergy per atom of the trimerized Al-ISAW with AFM ordering at several average interatomic 
distances d av . The nearest-neighbor interatomic distance between the atoms in the trimers 
and the bond angle formed by two adjacent trimers are optimized until the maximum forces 
on the atom are smaller than 82.4 pN. The total energy of the wire for N atom =2 given in 
Fig. |l| is also depicted in the column (B) of Table | for comparison. From Table H, one can 
see that before d av increases up to 5.2 a.u., where the ISAW in Fig. [I] has not changed to 
the linear wire yet, the ISAW forms a trimerized structure. We further comfirmed that the 
ISAW with AFM ordering is not a conductor but an insulator at d av =5.2 a.u. Thus, the 
Al-ISAW exhibits the Peierls transition just before forming a linear wire in a metallic state. 

In summary, we have studied the relation between the atomic and spin-electronic struc- 
tures of the Al-ISAW using FPMD simulations. We have observed that the Peierls transition 
occurs in the Al-ISAW with magnetic ordering: the Al-ISAW has a zigzag structure in a FM 
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state below d av =h.2 a.u. Just before the average interatomic distance increases to 5.2 a.u., 
the wire ruptures to form a trimarized structure with AFM ordering. The Al-ISAW changes 
from a conductor to an insulator upon breaking. The preference for a trimerized structure 
with AFM ordering in the distorted wire is attributed to the nature of its cr-symmetry bands. 

This research was partially supported by the Ministry of Education, Culture, Sports, Sci- 
ence and Technology, Grant-in-Aid for Young Scientists (B), 14750022, 2002. The numerical 
calculation was carried out by the computer facilities at the Institute for Solid State Physics 
at the University of Tokyo, and Okazaki National Institute. 
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